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STUDY OF A STEREO-ELECTRO-OPTICAL TRACKER SYSTEM 
FOR THE MEASUREMENT OF MODEL AEROELASTIC DEFORMATIONS 
AT THE NATIONAL TRANSONIC FACILITY 

Richard J . Hertel 
ITT Aerospace/Optical Division 


1 . 0 SUMMARY 

This study examined an electro-optical method to measure the 
aeroelastic deformations of wind tunnel models. The system studied is an 
electronic version of the photo grammetry method described by Brooks and 
Beamish. 1 In place of the film/ its processing and measurement/ an 
electronic camera system is used. Brooks and Beamish measured the 
recorded image positions of white target dots emplaced in the model. The 
electronic camera system does the same task in real time. It measures the 
coordinates of targets focused onto its photosensitive areas. 

This study shows the electronic camera system is capable of 
locating, measuring, and following the positions of 5 to 50 targets 
attached to the model at measuring rates up to 5000 targets per second. 

The targets must be either sources or efficient reflectors of light. 

This study analyzed, modeled, and measured the multi-target 
tracking performance of one of the two electronic cameras comprising the 
stereo pair. Although the conditions used for this study were those esti- 
mated to exist in the wind tunnel at the National Transonic Facility, NASA 
Langley, Langley, Virginia, it is emphasized tht a system of this type has 
broader application. This system can be used in other wind tunnels where 
the environmental and vibrational conditions are less severe than NTF. 

This system, can in fact be used for a rather wide range of three dimen- 
sional surface contour and displacement problems. 

This study considered the properties of the targets at the 
model, the camera optics, target illumination, number of targets, acquisi- 
tion time, target velocities and accelerations, tracker performance, and 
other factors affecting the measurement accuracy of target positions . 

Stereophotogrammetry has been tested at Langley Research 
Center as a means of measuring model deformations. Consideration of 
stereographic transformation was not included in this study inasmuch as 
the same or similar techniques used in the Langley tests will be applied 
when using the electro-optical system. 

The results of the study are presented here in terms which 
describe the eventual measuring system; its organization, installation, 
and operation. A baseline system design is given with comments regarding 

^Measurement of Model Aeroelastic Djeformations in the Wind Tunnel at 
Transonic Speeds Using Stereophotogrammetry, J.D. Brooks and J.K. Beamish, 
NASA Technical Paper 1010, Oct. 77. 
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alternatives. Some areas of technical uncertainty remain and these are 
noted. An analytical and exerimental section contains the information 
needed to support the baseline design and provides the basis for choosing 
among design alternatives. This report also describes the experiments and 
the demonstration, conducted at ITT-Aerospace/Optical Division, used to 
support the analytical work, to demonstrate feasibility of concepts, and 
to show working software and hardware. 


2 


2.0 


SYSTEM DESCRIPTION 


The system design described in this section is a conceptual 
one. Most features have only been treated analytically. Some features 
have been modeled in hardware, and in a few cases experiments have been 
performed to gather supporting data. 

The system described here assumes the use of small, self- 
luminous targets attached to the model. The targets are viewed by two 
electronic cameras located in the plenum of the wind tunnel behind optical 
ports in the test section wall. Operating modes provide for system 
calibration, target acquisition, target tracking and position measurement. 
Software dialog with an operator is used to select the system codes. 

Tunnel supervisory control is expected to identify when valid data logging 
conditions exist. This system then responds to an external command and 
records target positions while the model angle of attack is slowly varied. 
The recorded position information is analyzed later, off line, by the 
experimentor to extract the three dimensional target motions from the two 
camera stereo-optical data. Provision is made for a limited amount of 
quick look examination of the test data. Camera installation and siting 
are important considerations in the estimated system precision and 
accuracy. 

2. 1 Tunnel Installation 

* 

Figures 1 and 2 show the general arrangement of an aircrft 
model, electronic cameras, camera fields of view and coordinate systems 
used to relate three dimensional model coordinates to two dimensional 
image plane coordinates. 

The aircraft model is supported by the sting and is located 
along the center line of the tunnel. The model is movable in pitch angle, 
-11° to +19°, about a center of rotation located within the model. The 
model space coordinate system is designated (x" , y", z" ) and has its ori- 
gin at the center of rotation. The model support is made as rigid as 
possible but it does still bend and vibrate. The bending of the sting 
with a model attached gives rise to rigid body motions of + 0.5 degree 
amplitude in pitch, roll, and yaw about the center of rotation at frequen- 
cies of 8 to 50 Hz. Furthermore, the first bending mode vibrations of the 
model wing are expected to fall in the range of 25 mm peak-to-peak at 20 
Hz to 2.5 mm peak-to-peak at 200 Hz. 

The targets are sources of light, less than 500 pm diameter 
and flush with the surface of the model. The targets are either light 
emitting diodes (LEDs) or the polished end of a fiber optic light pipe. 

The electrical power input to each target is 100 mW for the LEDs. The 
fiber optic targets require greater power; depending upon fiber and 
coupling losses. The targets emit light in all directions, and at angles 
as great as 75° to the normal of the target. All targets must be within 
the fields of view of both cameras. These and ocher properties of targets 
having great importance to system operation are discussed further in 
Sections 3.2, 3.3, and 3.6. 
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There are two image dissector cameras used in this tracker 
system. Figures 1 and 2 show these cameras located 400 lun off the center 
line of the model, at positions upstream and downstream of the model. 
Camera locations for models having wing spans totally ins! ie the field of 
view are along the center line of the tunnel, amin, upstream and 
downstream of the model. In Figures 1 and 2, the emsetas lines off sight 
are aimed at a common point located at y" * -400 sun. The camera pointing 
, angles are (ct,f3,Y) * (125°, 90°, 35°) as shown in the figure. For the 

center line located cameras, the aim point is the canter of rotation. 

The system field of view is 800 by 1000 mm at the model for 
zero degree angle of attack and 1650 mm camera to model range. Figure 2 
illustrates an off center line camera installation . The field of view 
covers the fuselage and one wing of a 1.6 to 1.8 m wing span model. The 
boundaries of the field of view volume are set by the 34° x 34° angular 
field of view of the cameras and the 1650 + 450 mm depth of focus of the 
. camera optics. Greater depth of focus can be obtained by increasing the f 

number of the lens at the expense of increased target illumination power 
for the same data rates. The field of view can be changed by using a dif- 
ferent focal length camera lens. A shorter focal length lens gives a 
wider field of view and greater depth of focus at the expense of measure- 
ment precision at the model. A longer focal length lens has the opposite 
effect. Carried to an extreme, very short focal length lenses will have 
greater distortion while very long focal length lenses may have large f 
numbers. Section 3-1 on ortics discusses these considerations in greater 
detail. 

Each camera head, camera optics, and associated electronics 
are located inside temperature controlled pressure vessel. A vessel iso- 
lates the camera from the wide range of temperature and pressure con- 
ditions in the tunnel. The pressure vessel has a 10 mm thick optical 
quality window at one end. This window withstands the static tunnel 
pressure. A second optical quality window in the wall of the tunnel test 
section isolates the pressure vessel from the direct tunnel flow. A 
sketch of this arrangement is shown in Figure 3. 

The proper mounting of the cameras in the pressure vessel and 
the pressure vessel in the wind tunnel are important considerations * The 
mounting must minimize changes in the camera-to-model vector due to vibra- 
tion or temperature. The mounting must also allow for temperature control 
of the pressure vessel. 

Mounting point temperatures can range from -195°C to 95°C 
making both heating and cooling necessary. The camera and electronics 
present a 10-12 watt heat load. A satisfactory mounting is just as impor- 
tant as the proper target characteristics. While mounting details are 
outside the scope of this study, Section 3.6 discusses the precision and 
accuracy consequences of variations in the camera to model vector. 
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Figure 3. Camera Installation at the Test Section Wall 



2.2 


System Operation 


The system hardware is intended to operate upon the command of 
the system computer. Figure 4 shows a block diagram of the stereo-electro 
optical tracker system and the interface with the tracker system computer. 
Bach camera* uses a 16 bit, duplex, I/O channel, a 16 bit DMA ciannel and 
the associated command and control lines. LED targets are controlled via 
a serial I/O channel. If Camera Status is only a monitoring function, 
status data can be appended to the 10 bit video data. If command of the 
camera housekeeping proves necessary, another I/O channel will be needed; 
one channel for two cameras. 

The computer programs fall into three categories: 1) target 

acquisition, 2) target tracking, and 3) system calibration. Single 
camera examples of the first two types of software v/ere created during the 
course of this study and elements of all three types were used for the 
experiments and demonstration portions of this study. 

The flow chart in Figure 5 and the following sections explain 
the essential features of the system software. The main measurement 
program is TBCSS. TBOSS handles the operator dialog and the transitions 
between acquisition and tracking. Calibration i3 done another program 
using many portions of TBOSS* 

2.2.1 Target Acquisition 

The system software begins target acquisition by a call to the 
subroutine ACQDIR. This subroutine turns on all the LEDs and causes 
camera 1 to begin a raster scan of the field of view. Subroutine SEARCH 
moves camera 1 scan in x and y increments equal to one-half the image 
dissector aperture size. At each position in the field of view, the 
camera measures the scene flux and returns a numerical value to SEARCH. 
SEARCH calculates the difference between the signal and the neighboring 
background. If the signal exceeds the background by predetermined 
threshold value a target is found. The magnitude of the threshold value 
is based upon background level and signal-to-noise ratio needs of the 
system. (See Section 3.2.5 for a further discussion . ) 

Once a target signal is found, TGTVAL is called upon to vali- 
date the signal. TGTVAL checks the target list to determine if current 
camera 1 scan coordinates represent a new target or a repeat of an earlier 
discovered target. If no target is listed at or near the current coor- 
dinate, the signal-to-noise ratio above background is measured. If the 


*Appendix contains a paper describing the operating principles of image 
dissector tubes and cameras in tracker applications . The main text of 
this report assumes the reader has a basic understanding of this device . 
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Figure 4. System Electrical Block Diagram 




























Figure 5. Program TBOSS Flow Chart (Sheet 1 of 3) 
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signal-to-noise ratio meets the system needs, the current coordinates are 
added to the target list. If a target already exists at or near the 
current coordinate, or if the signal-to-noise is less than the minimum 
required, the target list is not changed and the acquisition scan is 
resumed. Acquisition scan continues until the next above threshold 
signal or until the entire field of view is scanned. 

The time, T, required to do a search is a function of the size 
of the field of view, the image dissector aperture, and the number of 
targets 

T » K , S + Nt 
1 1 2 

2i 

where d 

K-) is a constant relating the scan step size to the image 
dissector aperture size; typically 1<K-|<4. 

S is the size of the field of view at the camera faceplate; 

S = 10 x 10 mm. (See Section 3.1.) 

d is the size of the image dissector aperature; d = 200 x 
200 ym. (See Section 3.4.) 

ti is the time to access, measure and test the brightness 
value against threshold; ti ; 10“ 3 sec. 

N is the number of targets 

t 2 is the time to measure the signal-to-noise ratio and up- 
date the target list, t£ = 0.2 sec. 

T = (2 2 ) ( 10 2 ) ( 10" 3 ) + 50(0.2) 

( 0 . 2 2 ) 

= 10 + 10 seconds 

for 50 targets. 

Each camera can go through the acquisition sequence indepen- 
dent of the other camera; creating its own target list. The two target 
lists must, however, be reconciled. Both cameras, when measuring target 
positions must view the same target at the same time. Failure to measure 
moving or vibrating targets at the same time at each camera site leads to 
erroneous three dimensional position results • The target is at two dif- 
ferent three dimensional positions, neither of which has full stereo- 
graphic data . 

Two, reconciled target lists are obtained by causing camera 1 
to scan the field in acquisition as described and shown in the Figure 5 
flow chart. Before camera 2 begins its scan, all target LEDs but the 
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i th on the camera 1 target list are turned off. Camera 2, via SEARCH ant 
TGTVAL scans the field of view until the target is found or until the scan 
is complete. If no target is found, a message is printed. If a target is 
found it is positively known that target i, camera 1 corresponds to target 
j, camera 2. This procedure is repeated until all entries in camera 1 
target list have corresponding entries in the camera 2 target list. 


The average time to locate the illuminated target with camera 
2 , knowing only that it is from the camera 1 list, that it is the only one 
on and all targets are uniformly distributed over the field of view, is 


T 


1 


. K 2 S 2 t 2 

1 _I L 

2 d 


+ t 

2 


1 (2 2 ) go 2 ) (io~ 3 ) + 0m2 

2 ( 0 . 2 2 ) 


- 5.2 seconds. 


The time for fifty targets is 260 seconds. 

Knowing where to begin the camera 2 scan, at the top or bottom 
of its field of view will reduce this time by almost one-half. 

2.2.2 Target Tracking 

Following target acquisition, the system switches to the 
tracking mode for measurement of target positions. The computer sends the 
current target coordinates from the target list to each of the cameras. 

The tracker circuitry in each camera does the position measurement inde- 
pendently. Figure 6 shows the sequence of steps involved. In Figure 6, 
the target image shown as circular area, represents the optical point 
spread function seen in cross section. The image dissector aperture is 
square. Typical sizes are 35 ym full width half maximum for the target 
image and 200 x 200 ym for the aperature. Sections 3.1 and 3.3 provide 
reasons for these sizes. 

The tracker scan generator causes the image dissector camera 
to execute a rapid, four position track scan about the estimated target 
position. The scan amplitude in each of the two camera scan axes is typi- 
cally one-half the aperture size. (See Section 3.5.) The seven step 
sequence positions the center of the scan pattern to the center of the 
target.* The center of the scan pattern is eventually returned to the 
computer to update the target list of each camera and to be recorded as 
the target position. This sequence is repeated over and over for each 
target . 


♦Actually, the scan goes to a position having equal integrated brightness 
at each of the four sample positions. 
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Figure 6. 


Point Tracker Operating Principle. The X and Y 
Co-Ordinates are with Respect to the Scan Axes 
at the Camera 


There are several properties that describe the position 
measurement and tracking process. The first property is capture range. 

_1 Capture range is a vector measure of how far the target 
center can be from the center of the track scan and still 
permit the track scan to converge upon the target center 
position. Capture range is measured perpendicular to the 
camera line of sight. In first approximation, the capture 
range falls within the area covered by the convolution of 
the track scan with the target image. For the example 
described here , the capture range is shown in Figure 7 . * 

For successful operation, the target position estimate, 
after acquisition or when periodically returning to the 
target, must fall within the bounds of Figure 7. The size 
of the capture range sets a limit on the amplitude of the 
target vibrational motion.** For the Figure 7 case, the 
limit is approximately + 235 x.\J~2/2 pm or + 165 pm (+ 21 mm 
at the model ) . In some special cases , vibrations parallel 
to the camera scan axes, the limit is as large as + 220 pm 
( 27 mm at the model ) • 

2 A second descriptive property is the number of track scan 
cycles used to converge to the target center position. The 
number of cycles to converge is a function of the offset 
distance between target and the track scan. Experiment 
results of Section 3.5 show, for targets at the extreme of 
the capture range, 8 to 12 track scan cycles are needed. 

For target offsets less than 50 pm, one-fourth the aperture 
size, the convergence is completed in four or fewer cycles. 

_3 The third property is termed convergence criterion. 

Convergence criterion is the maximum acceptable uncertainty 
in the location of the target center. In the limit, for an 
image dissector tracker, the uncertainty in target position 
depends mainly upon the signal-to-noise ratio of the illu- 
mination falling on the photocathode. In practice other 
terms such as the deflection axis D/A converter resolution, 
target or camera vibration, seeing conditions through the 
boundary and shock layers also enter. The experiments in 
Section 3.5 show that + 0.5 pm is possible at the photo- 
cathode for electrical signal-to-noise ratios of 30:1 RMS. 
Vibration effects and seeing conditions are assessed in 
Section 3.6. 


* The capture range at the model depends upon the range to the model, the 
view angle, the attack angle and the position in the field of view. The 
size at the model can be estimated for zero angle of attack at the center 
of the field of view and 1650 mm range by multiplying the Figure 7 
dimensions by 125. 

**Section 3.3 shows target loss is expected to be infrequent. The system 
is programmed to skip to the next target on the list upon loss of signal. 
The dropout target is re-examined the next time through the target list. 
If all targets are lost the system switches to acquisition. 
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4, A fourth property used to describe the measurement and 
tracking process is the time required to do one cycle or 
more of the Figure 6 sequence. The tracker hardware is 
specifically designed to execute the tracking scan and 
signal processing algorithm. It causes the camera to scan 
the target/ it processes the video/ corrects the scan posi- 
tion/ and returns the corrected value to the computer. 

Figure 8 is a timing diagram of the tracker action for one of 
the two cameras. The camera goes to the +X scan position and integrates 
the signal for 20 ys. The time to store the +X value and go to the -X 
scan position is 5 ys. The video at -X is then integrated for 20 ys. The 
next 40 ys are used to update the estimate of the target X axis position 
and to send the camera to the +Y scan position. The signal integration 
sequence for + and -Y positions is the same as for the X position. Again/ 
40 ys are needed to update the estimate of the target X position and to 
send the camera back to the +X scan position. If the camera goes to the 
next target, 70 ys are needed to update the estimate of the current target 
Y position, store X, Y and send the camera to the estimated position of 
the new target. 

The time to do 1 cycle is ti « 130 ys. 

The time to do n cycles on one target is 


t(n) = (n) 130 ys + (n-1) 40 ys. 

The time to do N targets of n cycles each is 
T(N,n) » Nt(n) + (N-1) 70 ys 


= Nn 130 ys + N(n-1) 40 ys + (N-1) 70 ys 
For example if N = 50 targets, and 

n *» 3 cycles per target 
then t( 1 ) <=* 130 ys 

t (8 ) = 470 ys 
T(50,3) = 27 ms 

This is a rate of 1800 targets per second. The highest rate 
would be for n = 1 or 5000 targets per second. 


Not all targets can be measured to an uncertainty within the 
convergence criterion in a single cycle. Neither do all targets require 
3, 8, or even 12 cycles. It is advantageous to make the tracker adaptive. 
The tracker remains at a target until the convergence criterion is met or 
until a maximum number of track scan cycles have been completed; then goes 
to next target.* 


*This action was included in the demonstration software. 
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2 . 2,3 System Calibration - The system user needs to know the trans- 

formation function between model space coordinates and tracker system out- 
put data. More specifically the user needs to know how the system differs 
from the mathematical model built into the data analysis software. 

The transformation program currently being used at NASA 
Langley for photogrammetry corrects for symmetric and anti -symmetric lens 
and camera distortion in a manner similar to the one used to determine 
viewing angles and camera distances. That is, it is a least squares coef- 
ficient fitting technique using known position of reference points. It 
is thought that some reduction in computational effort could be achieved 
by measuring the mapping function of the camera, lens, pressure vessel 
window and perhaps even a test section window in calibration fixture out- 
side the tunnel. Section 3.6.1 describes an interpolating distortion 
correction method that i3 estimated to provide correction to + 0.5 yrn at 
the camera phot&cathode. The camera, lens and windows with a correction 
algorthium and data base would then yield a system having a linear mapping 
function + 0.5 pm X magnification between the object plane and camera data 
output. 


The data base is obtained by mapping the field of view of each 
of the two cameras. It is assumed the distortion is a slowly changing 
function of target position within the field of view and that a finite 
number of points are sufficient to describe the function. (Section 3.6.1 
shows that 625 points are needed.) 

The test fixture must have means to move a point source of 
light about the field of view. Such a means is an array of LEDs located 
at known positions on a flat plate. There must be provisions to assure 
the plate is 1650 mm from the camera and that it is perpendicular to the 
optical axis. Variations on the acquisition and tracking software gather 
the distortion map data. An analysis program is used to create the 
calibration data base. 

The camera and lens distortion calibration map is either pro- 
vided with the recorded tunnel test data or could be used to correct the 
recorded data before the data is given to the experimenter. Since the 
tracker system must operate in real time when recording data, it is better 
to do distortion corrections off line, after all data is collected. 

Estimates of the optical window and the optical path seeing 
conditions upon the measurement accuracy are given in Section 3.6. 

It is still necessary to calibrate mm at the model as a 
numerical value in the software. The stereographic method being con- 
sidered here is a non-metric method in the projective geometry sense- 
This means that sizes or scale factors are not preserved and thus must be 
calibrated. It also means that the exact viewing angles of both cameras 
and their separation distance need not be known. The method 1 does need at 
least 6 and preferably as many as 20 reference points surrounding and 
throughout the volume of the combined two camera field of view. The 
reference points must not be in the same plane. 
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This portion of the calibration must be done in the tunnel 
with camera installed. The 0.6 to 20 reference points can be specific 
targets of known location on the model. This calibrat on should be done 
at test temperature and pressure as part of the data logging process. 

2 . 3 System Computer 

The computer used in this system must be capable of supporting 
the operation of the two image dissector cameras , the target LEDs, a 
magnetic tape recorder, a CRT operator terminal, and printer. It may also 
be desirable to have a card reader and a floppy disc unit for ease of 
programming and program modifications. 

During acquisition mode, each of the two image dissector 
cameras are controlled via 16 bit, duplex interfaces. Instruction words 
are used to set internal camera conditions and to select or not to select 
the tracker mode. LEDS operate via serial inteface. Instructions are 
sent only when a change in operation is required. 

Two words are required for each (X,Y) coordinate. If one 
coordinate is unchanged, only the one word is needed. The rate during 
acquisition is 1000 coordinates per second, set mainly by software execu- 
tion time. 


Scene flux measured at selected coordinates is returned to the 
computer as single 16 bit words. The rate is again 1000 measurements per 
second. 


During track 'and measurement mode, the exchange between com- 
puter, LEDs and the two cameras occurs at a 40 K word rate. The computer 
must, at the same time, service the magnetic tape drive at an average 20 K 
word rate, monitor the "log data" command from the outside, update the 
target list, and coordinate the actions of two cameras. The magnetic tape 
unit operates via a direct memory access interface. The tracker circuitry 
has coordinate storage buffers containing the target list. These buffers 
communicate with the computer via direct memory access . Thus the CPU 
should have ample time to observe and direct the system activities without 
being burdened with all the data transfers. 

During calibration, the computer ^ammunicates with the cameras 
as already described. The target light sources are now those on the 
calibration fixtures. 

The acquisition software is written in a high level language 
such as Fortran. The driver modules for the cameras and light sources are 
written in assembly code. The track subroutine is in Fortran. The track 
scan and video processing algorithm (Figure 6 sequence) is in hardware at 
each camera . 


No estimate has been made of memory requirements. The experi- 
ments and demonstration were conducted under a real time operating eystem 
and Fortran on a computer having 24K words of memory. The main tracking 
program plus operating system used 16K words. For the experiments, the 
track scan and video processing algorithm was in Fortran. 




3.0 ANALYSIS AND EXPERIMENTS 

The design of a stereo electro optical tracker system used to 
measure aeroelastic deformations of wind tunnel models must consider the 
following essential elements. 

1 . The characteristics of the targets attached to the model • 

2. The illumination of the targets and the model. 

3. The optical properties of the media between the model and 
the electronic cameras. 

4. The characteristics of the electronic cameras. 

5. The effects of the tunnel environment and system installa- 
tion upon operation. 

3 . 1 Optics 

The field of view i3 800 x 1000mm at 1650mm camera to model 
distance and zero degree angle of attack. Image dissector camera tubes 
are available in three sizes having photocathodes of 17, 25, and 44mm 
diameter. Figure 9 shows these photocathodes, the image format size, and 
appropriate focal length lens. 

The model can be moved +/- 450mm about the nominal 1650mm 
distance. It is necessary to maintain focus over this range of distances. 
Operation off focus increases the radius of the point spread function of 
the lens at the image plane. The change in radius due to off focus opera- 
tion can be estimated by 

AD. 

r » — r— tan (a) 

where 

a =* sin - ^ ( 1/2f ) 
f «* the f number of the lens. 

. Large f numbers provide good depth of focus but are very inef- 

ficient collectors of light. A compromise must be achieved between depth 
of focus and lens aperture. 

The magnitude of ADj., the change in back focal plane position, 
depends upon the focal length of the lens. 

(FL) 2 


for FL/Do « 1 


; 

\ 

r 
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IMAGE 10 x 10 nun 15 x 15 mm 25 x 25 mm 

FORMAT 


LENS 

FOCAL 16 25 41 mm 

LENGTH 


AD . 

i 

FOR AD =900 .091 mm .22 mm 0.60 mm 

o 


Figure 9* Camera Image Formats and Lens Focal Lengths for Three Sizes 
of Image Dissectors. AD. is the change in back focal plane 
position caused by +450 mm change in object position for an 
object at 1650 mm. 




Calculated values for Dj_ are included in the table in Figure 8 


If A r is scaled to the size of the field of view, that is 

FOV 



Ar 

and FL 

then the required f# is 

The results favor the use of a sma.V, format tube and short focal length 
lens . 


The light collection efficiency of an f/7.3 lens is only 15 
percent of the collection efficiency of the f/2.8. It is shown in Section 
3.2 that such a loss is not desirable. 

Qg^The growth in the point spread radius for the f/2.8, 16mm lens 
is r = tan (10.3) = 8fjm. This value is tolerable. The electron 

optics of the image dissector camera tube exhibit a 5- 10pm radius point 
spread. T.n f/2.8 lens has a finite point spread radius also, about lym 
theoretically but more likely 5-10 ym in commercially available unit of 
conventional design. Thus the 8 y m change in point spread radius should 
be compared to the 10-20ym radius point spread. It is shown in Section 
3.5 that absolute size of the point spread is of lessor importance than 
changes in size. 

The important properties of the tracker lens are low distor- 
tion, 0.5%, flat focal plane, + 10 ym symmetrical abberations, low scattered 
light, and mechanical stability. 

3 . 2 Targets 

Targets cannot be considered without also considering some 
aspects of target illumination and tracking method. The model can move in 
pitch, roll and yaw over a substantial fraction of the tracker system 
field of view. Thus the target appearance can change substantially. The 
ideal solution is a combination of targets, illumination, and tracking 
method which yields a position measurement accuracy independent of changes 
in 

1 . Target size - magnification 

2. Target orientation - rotation 

3. Target shape - distortion 

4. Target illumination - magnitude, uniformity 

5 . Target contrast - stray light 

6 • Target focus - depth of focus 

7 . Target viewing angle - off normal view 

The target and its illumination must also be optically efficient in terms 
of power input to the illuminator for the signal needed at the camera. 
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3.2.1 Point versus finite sized targets . - There are two general 

classes of targets? point and finite. Point targets are not resolved as a 
disc image by the camera and optics. Finite targets have a definable 
size, shape and even structure. 

Point targets can be defined quantitatively as having an image 
whose size is less than one- third the size of the camera optical spread 
function. The section on camera optics estimated the radius of the spread 
function to be 10-20 ym plus a 0 to 8 ym term for depth of focus effects. 
Using the dimensions of Figure 1 and a 16mm FL lens, the smallest demagni- 
fication is about 75 at a range of 1200mm. The size of the target, at the 
model, is therefore 2 x 750 ym/3 = 500 ym or less diameter. 

Target position for point targets is measured by making four 
samples in the vicinity of the target image as shown in Figure 6. The 
pointion of the four samples is adjusted until all have equal signal 
amplitude. 


Numerous types of finite targets can be considered. In 
general finite targets require more image samples than point targets 
because of the range of sizes and possible orientations. 

Two simple examples of finite targets are the white circle in 
a black surround or the black and white wedge. 



Figure 10. Two Examples of Finite Targets 

Since the model can move in range and pitch angle, these 
targets in the extreme can look like those shown in Figure 11. A scanning 
method which is independent of these size and shape changes is a raster 
whose dimensions are greater than the maximum possible target size. The 
scanned data is then used to compute the center of brightness ; assuming a 
uniform target. This method is time consuming if done by computer or 
complicated electronically if done by a dedicated hardware processor. 
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Figure 1 1 • Finite Targets Viewed at Large Angles 

A less complex or time consuming method uses a circular scan 
pattern. The method seeks the black white edge crossings and adjusts the 
center of the scan circle until two pairs of edge crossings are each 
180° apart. This method is sensitive to target rotation. The rotation 
must be known or measured before proper adjustment of the scan center 
occurs . 


In general, finite targets require a tracking an measurement 
method insensitive to size and shape changes. Finite targets are best 
used where some useful information is conveyed by their size, shape, or 
pattern. 


3.2.2 Active versus passive targets . - Targets can be either active 

sources of light or passive reflectors of light. Active point targets 
such as light emitting diodes are optically very energy efficient. All 
the energy contributes to the source intensity. Fiber optic light pipes 
are a bit less efficient due to losses along the fibers and at the input 
coupling. Nonetheless, optical fibers viewed end on are potentially good 
targets . 

Passive point targets such as diffuse painted patterns, are 
optically very inefficient in this installation. Only a small fraction of 
the incident energy is reflected to the cameras. Passive point targets 
can be illuminated by a spot light that moves in step with the tracker as 
the tracker moves from target to target. Such a scheme would require 
either premapping of all nominal target positions versus angle of attack, 
or, real time computation of illuminator pointing angles. It is still 
necessary to illuminate an area larger than the passive target since an 
allowance for pointing errors must be made. 


The flux onto a target, at normal incidence is given by 

4 it f 2 SNR 2 M 2 Ke£ £ 

F “ V S At 12 
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A = lx 10“ 6 m 2 , tlie a ^ea of a point target 

f = 2.8, the focal ratio of the lens 

SNR = 30 , the signal-to-noise ratio needed according to 

Section 3.5.2 

M = 125, the demagnitifcation between target and image 

K = 2.3, the noise factor of the camera 

e = 1.602 x 10”^® coulombs/electron 

e l e 2 = re fl ect i° n and lens transmission factors 
S =* 200 pA/lumen, the camera sensitivity 

At =» 20 ys, the time to make one brightness sample 

) F = B A = 4 ^(2.8) 2 (30) 2 (125) 2 (2.3) Cl. 6 X 10~ 19 ) (1) (1) 

° (200 x 10~ 6 ) (20 X icf 6 ) 

= 0.13 lumen 



The target is only 1 x 10“® m 2 in area, yet the flux is 0.13 lumens. The 
illuminance must therefore be 1.3 x 10^ lumen/m 2 . This is direct 
sunlight. For flood illumination of the 1m 2 field of view at the model, 
this is at least 8 Kw of quartz-iodide lamps in housings with very effi- 
cient reflectors. For a moving spot of illumination using a 10 
cm 2 area this is about 250 watt of quantz-iodide illumination with an f/2 
projector lens. 


A corner cube reflector used as a passive point target' can be 
energy efficient. A collimated beam of light to a corner cube is 
reflected back along the same path to the source. A corner cube target 
does not exhibit the depth of focus effects previously mentioned because 
the source appears to be located at infinity. There are two major 
problems with corner cubes; size and view angle. It is not clear if 
millimeter sized cubes are available at an economical cost. Corner cubes 
do not work at large angles ( 45°) of incidence. 


In general , 


active, point targets such as LEDs are preferred. 



27 


3*2.3 Target contrast . - Contrast between target and surrounding 

model is an important property. A point target tracker can, in principle, 
operate on low contrast targets. In practice, the background signal 
should be eliminated to minimize the dynamic range of the electronics and 
to remove the effects of background flux changes. 


Let F 0 , F-j , and F 2 be the average values for the background, 
sample 1, and sample 2 signals for the Figure 6 tracking method. 
Associated with each sample is a noise n Q , n^, and n 2 , each proportional 
to the square root of their respective signals. The track scan position 
error signal, E, is: 


= C, 


< F 1 + V - < F 2 + F o ) 
(F 1 + V + (P 2 + F o } 


“ C, 


F 1 - F 2 


p i + F 2 + 2 F q 


It is possible to remove the 2F 0 term by making a measurement 
of the background signal. This additional measurement will add 20 to 60 
ys to the single cycle time shown in Figure 8. 

Removal of the 2F 0 term does not completely solve the 
background problem. At null position, F-j = F 2 , the finite noise asso- 
ciated with each of the samples becomes important. The noise, e, about 
the null position is equal to 

(2n 2 + n 2 + n 2 )* 5 

e ~ C 0 12 

1 Pi + f 2 + <an| + + „!)■> 

If F 0 = F-j, 50 percent contrast, then the position noise about null is at 
least V 2 " times greater than with zero background. As F 0 becomes larger 
than Fi , a still lower contrast, then nQ, the noise in the background 
signal becomes the dominant* term. 

If F 0 <_ 0.1 F^, a condition corresponding to >90 percent 
contrast, background measurement and subtraction is not necessary. If 
F 0 = 0.1 F'j , the noise about null is only 5 percent greater than the zero 
background case; an acceptable condition. 

Tracker control loop gain depends upon the 2F 0 term in the 
error signal equation. For 0 <F o £0.1 F^j , the change in loop gain is less 
than 1 dB, an acceptable value. (See Section 3.5.) 

The limits on absolute background brightness and the needed 
target intensity can be estimated for a system using LED targets . Section 
3 . 5 shows that the required signal-to-noise ratio for the difference bet- 
ween target and background is 30:1 RMS. This signal-to-noise ratio allows 
the center of the target to be located at + 0.5 ym at the model for zero 
angle of attack and center of the field of view. 


*Image dissectors will respond to the 120 Hz ac flicker from incandescent or 
fluorescent lights. This can increase the b noise term considerably beyond 
the shot noise value. 
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Let I 0 and Ii be the photocurrents entering the aperture of 
the image dissector due to the background and the target plus background. 
It is required that the signal-to-noise ratio of Ii“I 0 be greater than or 
equal to 30. Since the image dissector is a shot noise limited device, 
the noise in each of the two photo currents is 


where 


*1 

e 

k 

At 


ekl 
c 

. At 
eklj 

At 


amperes 


amperes 


1.602 x 10“ Coulombs/electron 
2.5, image dissector noise factor 
sample time, 20 ys in this system. 


SNR = 


I, - I 
1 o 


I- - I 

1 o _ 

(i 2 + i 2 ) 1 * ek 3 (I + I )** 
° 1 At 1 ° 


These equations assume no ac modulation of the background. If there is ac 
modulation, i Q becomes the rms value of the modulation. As calculated 
here the equation gives the greatest dc I Q value that is tolerable. The 
noises are assumed to be uncorrelated and to add in quadrature. 


Figure 12 is a plot of I-j-I 0 , the required target photo 
current versus background photo current I 0 for SNR = 30 and At = 20 ys. 

An additional abscissa scale on the plot shows the maximum background 
luminance at the model for 2854°K tungsten illumination, f/2.8 lens and a 
200 x 200 ym aperture in the image dissector. General office illumination 
is of the order of 100 lumen m”2 steradian“1. 


The extra scale on the ordinate axis gives the target source 
intensity to cause the signal photo current. This intensity is based upon 
an f/2.8 lens and a red LED. Red LEDs produce 1 x 10“^ watts/steradian at 
60 mA, current, 25°C, and normal incidence. 

Undoubtedly, there will be other sources of light used within 
the tunnel. These sources could be incandescent, fluorescent, or laser. 
They may even be turned ON and OFF for various reasons and will probably 
operate from the ac line . It is prudent to make the tracking system as 
insensitive as possible to these sources. 

There are several means available to achieve this end. First, 
make the targets the most intense source within the tunnel . LEDs can be 
pulsed to intensity levels 20-100 times greater than their dc, 25°C 
ratings for 1 percent duty cycle at 300 Hz pulse rate. In this system, 
the target diode would be turned ON only when the particular target posi- 
tion is being measured. This will require a small amount of circuitry in 
the model in addition to LED controller described in Section 2.2. 
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Further reduction in system susceptibility to background illu- 
mination can be obtained by using narrowband optical filters on the 
cameras. If a filter with a passband of 50 nm centered on the 655 nm LED 
emission wavelength is used, the insertion loss will be about -3 dB. The 
out-of-band transmission is usually less than -60 dB. 

With this kind of optical filter, background light due to 
tungsten sources would be reduced by -13 dB, due to cool white fluorescent 
sources by -19 dB, but only 0 to -3 dB against HeNe laser sources. 

Added discrimination against HeNe laser light can be obtained 
by using a narrower bandpass and accepting a greater insertion loss or 
using the less efficient green LEDs. The problem is the emission wave- 
lengths of red LEDs, and particularly the high efficiency red LEDs, 
overlap the HeNe laser wavelength. 



Figure 13. LED Relative Intensity Versus Wavelength 

If the narrow band optical filter does not provide sufficient 
discrimination between target and background, a phase coherent, modulated 
light system could be considered. Each of the target LEDs would be modu- 
lated at a frequency in the range of 5 MHz. Since the image dissector 
tube will pass this frequency, a tuned amplifier, phase sensitive detector 
would replace the conventional preamp. The system would accept optical 
radiation within the optical filter passband having an amplitude modula- 
tion of 5 MHz and the proper phase relationship to a reference signal . 

This is a complicated approach and is an extreme but effective measure . 
This approach avoids the 20 to 60 ps time penalty associated with 
measuring the background. The 50 percent duty cycle of the 5 MHz modula- 
tion doubles the needed target intensity. 
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3.2.4 Angular distribution of the target light . - An almost 

hemispherical distribution of light from the target is essential. 

For a flat plate at zero angle of attack the cameras view 
targets at a 35° angle from normal incidence. As the flat plate changes 
angle of attack, the 35° angle sweeps over a range of angles from 16° to 
54°. The camera has a rather wide field of view, 34°, therefore a nominal 
35° viewing angle at the center of the field of view becomes 18° to 
52° at the sides of the field of view for zero angle of attack and -1° to 
71° for 19° angle of attack. The greatest viewing angle is about 72° for 
19° angle of attack and occurs in the corner of the field of view. If 
targets are mounted on curved surfaces, this can lead to still larger 
viewing angles . 

For a target having a Lambertian distribution, the intensity 
varies as cos0, where 0 is the viewing angle from the normal. At 0 * 70 
the target intensity is only one-third of the normal value. LEDs are not 
Lambertian. Intensity at 70° is about 10 to 20 percent of on axis inten- 
sity. Stated differently, the brightness required for off axis viewing 
will determine the target source intensity. This fact alone will likely 
dictate pulsed operation of the LED targets . 

3 . 3 Target Dynamics 

Vibrations of the model will occur during tests . The vibra- 
tions will primarily be rigid body motion of the model on the sting and 
bending motion of the wings. Targets attached to the model will therefore 
vibrate about some average position. 

Table 1 gives the amplitudes and frequencies of the largest or 
most rapid vibrations as provided by NASA for this study. The motions at 
the model were transformed to motions at the camera image plane using the 
parameters 


*o 

= 946 

mm 

a = 

125° 

Y o 

= - 400 

mm 

3 = 

90° 

z o 

= -1352 

mm 

Y = 

35° 


and a back focal distance of 16.157 mm. These are the conditions 
illustrated in Figures 1 and 2. 

Other motions of the targets , such as those which occur with 
angle of attack changes, take place much less rapidly or are of very much 
smaller vibrational amplitude. 

All the motions given in Table 1 are excursions about an 
average position. The largest distance is 2 x 63 pm = 125 pm peak-to-peak 
at the camera for the 8-50 Hz yaw motion. The capture range of the 
tracking system should be larger than this dimension to avoid the need to 
periodically reacquire a lost target. Reacquisition is time consuming in 
software. It adds to the complexity of the tracker circuitry if done in 
hardware. For targets at the edge of the field of view and at 19° angle 
of attack, the vibration amplitude at the camera can be 25 percent larger, 
or 155 pm. 
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Table 1. Target Vibrations at the Model and at the Camera 


Rigid Body +0.5°, 8 to 50 Hz at 0° Angle of Attack 
I Roll | Pitch Yaw 


Target At 


-508 

-610 

+2.54 


-508 + 5.3 
-610 
0 


Units 


Motion at 
Camera 


Wing Bending, Target at x" = -508, y" -610, z" = 0 


Motion at 
Camera 



Az" = + 12.7 mm 
20 Hz 

Az" = + 1.3 mm 
200 Hz 


Units 


• 



+49 
+ 19 

+4.9 
+ 1.9 


y m 
y m 



qy. " tut ^pgp^ssw'' — 


The capture range shovm in Figure 7, for a 200 x 200 m aperture 
in the image dissector camera tube, has an assured capture circle of 

(435 - 200) x iC£ « 165 pm radius, 

2 

Thus no need for reacquisition software or hardware is expected. 

Should the vibration amplitudes be greater than those presently 
assumed, the U3e of a larger aperture will solve the problem at the expense of 
less tolerance to background light and a greater number of cycles to converge 
for large offsets . 

3.4 Image Dissector Camera 

The cameras used in this system are modifications of a commer- 
cially available image dissector camera and digital interface unit. The hard- 
ware tracker is an additional item of electronics, A description of image 
dissector tubes as trackers is in the Appendix. 

3.4.1 Imrja Dissector Camera . - The camera head contains the lens, the 
image dissector tube, the shielded deflection and focus assembly, voltage 
divider, .and ?,rt auplifxer . The deflection yoke driver, focus current regula- 
tor, nigh ' v qi power supply and tracker hardware are contained in an 
electro... o uni x located near the camera head and inside the pressure isola- 
tion vessel shown in Figure 3 • Other items of housekeeping circuitry such as 
the temperature controller and the digital to analog converters are also 
located with the camera in the pressure isolation vessel. 

The tube used in the camera has a 200 x 200 pm sampling aperture. 
The photocathode is a trialkali type, S-20 on a 7056 glass window. 

The camera operates static focus and has 15 bit addressing in both 
x and y axes. The addressable field of view is 12.5 x 12.5mm. One LSB in the 
deflection is 0.4 pm at the faceplate. 

The camera requires +/-15 volts supplied by the digital interface 
unit located outside the wind tunnel. 

3.4.2 Digital interface unit . - All communications between a camera and 
the computer take place through this unit. It decodes the computer commands 
into operating modes of the camera, tracker, and housekeeping circuitry. 

An auxiliary circuit located in the model provides the on/off 
control of the target LEDs. This circuitry occupies about 500 cm 3 and 
requires 7-10 watts, most used to light the LEDs. 
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3.5 


Tracker Characteristics 


During this study/ a single camera/ multi-target/ point tracker 
system was assembled using existing laboratory test equipment and hardware. 

The camera was a general purpose image dissector camera and its digital inter- 
face unit. A computer was programmed to run the camera. The track scan and 
video signal processing algorithm was written in software. This system 
contained many of the features that will eventually be incorporated into a 
stereo electro-optical tracker system. Section 3.8 describes the demonstration 
of this system that was made at ITT-Aerospace/Optical Division/ Fort Wayne/ 
Indiana. This section presents results of various measurements made on the 
system? measurements made to quantify the performance. 


3.5.1 Error detector . - The camera plus tracker logic is a feedback 

system. The system begins with a target position estimate/ measures the error 
in the estimate/ and corrects the position estimate. A key element is the 
position error detector. In this system it is not one discrete piece of hard- 
ware. but rather a term describing the results of scanning the target image 
and extracting a measurement of its position relative to the track scan. 

Figure 14 shows the measured value for the normalized error detector function 


E(x) 


F. (X) - F, (x) 
F~ (x) + F ^ (x) 


y * const 


This function relates the measured imbalance in the target brightness samples 
to a position difference between the center of the track scan and the center 
of the target. When this distance is sufficiently small, the target position 
is said to be measured. This curve can also be used to relate the video 
signal-to-noise ratio to an equivalent position noise. The curve is linear 
for small position offsets near null. Away from null the slope of the curve 
decreases. For distances greater than 50 ym from null and up to the limits of 
the capture range, the slope is near zero. 

The slope of the E(x) curve is shown in Figure 15 and is con- 
sidered to represent the sensitivity of the error detector. A change in error 
detector sensitivity is also a change in loop gain of the feedback system. 
Changes in loop gain carry implications for system transient response and sta- 
bility. For the purposes of these and other experiments, the loop gain was 
adjusted for best settling and minimum over shoot near null. This gain 
adjustment was done in software. 

The non-linearity of the E(x) function for large x is a nuisance. 
It slows the convergence to null for large offsets as will be seen later. 

More troublesome are phenomena which cause the slope near null to change. To 
assure system stability, the maximum loop gain must correspond to the steepest 
slope through null. Too high gain or an increase in slope could lead to 
oscillations about null. After the loop gain is set, subsequent loss of error 
detector sensitivity leads to increased convergence time. 
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The shape of the error detector curve and consequently its slope 
near null is affected by changes in the target image. A larger image will 
cause a reduction in slope. Magnification changes between the model and the 
camera or depth of focus effects can cause changes in the size of the target 
image. The design approach has been to make the target image smaller than the 
optical point spread function and thereby minimizing the magnification 
effects. The depth of focus effects have been minimized by use of a f/2.8 
lens, see Figure 16. A slower lens gives a greater depth of focus but at the 
expense of target intensity. There is some advantage to creating a larger 
point spread function. A larger point spread will tend to linearize the E(x) 
function over a greater range and thereby reduce the number of cycles needed 
to converge to null. However, such action requires the optical signal-to- 
noise ratio be increased in order to hold the same noise equivalent position 
error . 


The targets form an important part of the system design and have 
an effect upon performance. The experimenter and model builder will need to 
consider this fact. 

3.5.2 Noise equivalent position . - Noise equivalent position is a fun- 
damental property of the system relating to measurement precision. It is a 
way of expressing system noise sources in terms of a single position uncer- 
tainty at t'he faceplate of the camera or at the model. Some of the important 
noise sources are the shot noise in signal, the D/A and A/D quantization 
noises, and camera and lens distortions. The effects of vibration and seeing 
conditions are treated separately under Section 3.6 on accuracy. Camera and 
lens distortion are presumed to be fixed and measurable. 

Figure 17 is a graph of position uncertainty as a function of 
signal-to-noise ratio of the target signal. The target signal includes both 
shot noise and A/D quantization noise. The position uncertainty is asymptotic 
to + 1/2LSB* or + 0.25 ym. A signal-to-noise ratio of 30:1 yields a + 0.5ym 
position uncertainty. This curve resulted from measured performance on the 
hardware and software used for the demonstration. 

3.5.3 Measurement cycles for null convergence . - As mentioned in 3.5.1, 
the error detector function is non-linear away from null. This implies that 
several track scan measurement cycles are needed. Each cycle contributes, 

an increment of correction, as the error is nulled to the noise limit. 

Figure 18 is a plot of the measured number of cycles needed to 
converge +1LSB or + 0.5 ym of the null position versus the distance from null. 
The 0 dB gain condition represents no change from the loop gain that provides 
the fastest settling and minimal overshoot about null. The other curves are 
deliberate reductions in the loop gain to illustrate the effect upon con- 
vergence time. The time to do n cycles is 

t n = (n) 130 ys + (n-1) 40 ys 
based upon the Figure 8 timing diagram. 


*The demonstration system provided this data and used 13 bit deflection over 
4.1 mm or 0.5 ym/count. The system described in Section 3.4.1 is the one 
to be used in the tunnel installation. It will use 15 bit deflection over 
12.5 mm or 0.4 ym/count* The signal-to-noise ratio dominates the + 0.5ym 
position uncertainty. 


dE(x) yrrf 


I 




^mj C pAr ITY 0F VIP 

mAL page j S j> ( ,ur 


•FOCUS 


125]im 


250ym 


■375ym 


500ym 


■50.00 


-25.00 


25.00 


50. 0C 


MICROMETERS 


Figure 16. 


°lope Changes in the Normalized Error Detector Curve 
Changes in Focal Plane Position. For an f/2.8 Lens, 
1650 ± 450 mm Change in Object Distance Imolys A ± 4 
Change in Focal Plane Position. 




-l.ldB 


01 S310A0 


Figure 18. Cycles to Converge to Null Versus Offset from Null 
and Loop Gain (1 count = 0.5 ]jm) 



The stepped, even sawtooth appearance is due to a non-linear 
change in loop gain implemented in the software tracker used for these tests. 

When the position correction is calculated to be greater than 20 counts or 10 
pm, the correction actually made is 40 counts or 20pm. This scheme initially 
overshoots but cuts the number of cycles in half. 

Since a target can be expected anywhere within the capture range 
when the tracker returns to re-examine its position, quite a few cycles could 
be required to reconverge some targets. Tracker hardware design should allow 
a degree of adaptation. It should require the target measurement cycles to 
continue until the correction value is less than or equal to +1LSB or until 
some maximum number of cycles have been executed. This will assure full measure- 
ment precision with minimum time expenditure, plus a means to identify less 
accurate measurements , yet not be trapped forever while attempting to converge 
high amplitude, rapidly vibrating targets. 

3.6 System Precision and Accuracy 

This section discusses the precision and accuracy of the system 
concept. Precision is taken to mean the ability to determine small changes in 
a target position using a stable, repeatable scale of measure; a scale that is 
consistent throughout the field of view. The accuracy of the system is taken 
to mean the degree of agreement between the system scale of measure and an 
external absolute scale such as millimeters at the model. If good system 
design yields a sufficiently precise system, then a suitable calibration [ 

method and data reduction should give a commensurate accuracy. * 

The limitations upon the precision of the system will be covered 
by examining the contributions of the several parts. The major divisions are 
the camera and its lens , the dynamical motions of the targets , the dynamical 
motions of the camera installation, and the seeing conditions. 

Ideally the results should be described in terras of dimensions at 
the model. However, there is no single factor that relates incremental posi- 
tion changes at the model to image changes at the camera. The relationship 
depends upon the location of the target in the field of view, the viewing 
angle of the camera, the focal length of the lens, the range between model and 
camera, the angle of attack of the model, and vector direction of the position 
change. 


For estimation purposes the model can be considered a flat plate 
at zero degrees angle of attack, 1650mm range and viewed using a 16mm focal 
length lens. For a target located at or near the center of the field of view 
and viewing angles (agy) = (125,90,35) (See Figure 1), the model to image rela- 
tionship can be described by 3 de-magnification numbers, Mx" My" Mz" . 

MODEL CAMERA 


SPACE 


SPACE 

FACTOR 

Ax" = + 1mm 

Ax 

= + 8.0um 

Mx" = 125 

Ay" + 1mm 

Ay 

+ 9 . 8ym 

My" = 102 
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Az 

+ 5.6ym 
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3.6.1 Camera and camera lens effects upon system precision . - The camera 

and lens act in several ways to limit the system precision. The simplest 
limitation to understand is the one due to small but finite sized deflection 
increments. The demonstration system, described in Section 3.8, used 13 bit 
deflection D/A converters over a 4.1 x 4.1 mm field of view. The minimum step 
size was 0.5 pm at the camera faceplate. The system design concept presented 
in this report uses a 10 x 10 mm working field of view. Fourteen and 15 bit 
converters provide 0.8 to 0.4 Pm increments respectively while covering a 12.5 
x 12.5 mm scanable field of view. Such converters are available for the 
camera . 


D/A converters and cameras have temperature dependencies which 
affect the size of the area scanned and the stability of the scan center 
within the field of view. A change in scan size is equivalent to a change in 
scale factor at the camera faceplate. The magnitude of the error is equal to 
the distance from the center of the field of view times a +20 ppm/°C tem- 
perature coefficient. Values range from zero at the center to +0.12 pm/°C at 
the edge of the field of view. 

The centering term is +10 ppm/°C of the full scale size or +10 x 
10" 6 x 12.5 mm = +0.12 Um/°C. All coordinates within the field of view are 
equally affected by this factor. The pressure isolation container used to 
house the cameras must have temperature control because of the -195°C to 95°C 
tunnel environment. The temperature coefficients require the control of the 
temperature, at certain key locations in the container, be the order of plus 
minus one to two degree Celsius. 

The mechanical packaging of the camera and its lens must prevent 
shifts of the tube relative to the focal plane or coil system. The need is 
for less than +0.5 pm instabilities, a value consistent with the needs of 
precision star trackers used for space flight. 

The camera and lens both introduce distortions of the field of 
view. Figure 19 shows the magnified distortion of an uncorrected image 
dissector camera without lens. It is assumed the lens distortion is com- 
parable or less in magnitude. The dotted boxes are on a 1 mm grid and cen- 
tered on the true position . The distance from the center to the side of the 
smaller boxes is + 25 pm; + 40 pm for the larger boxes. 

This distortion is in principle correctable using a lookup table 
and interpolating between table entries. A key question is how large must the 
table be in order to correct + 0.5pm? 

Consider a group of four boxes in Figure 19 whose centers are 
denoted by A, B, C, D. The distance between AB is 2 mm and the same holds for 
CD. Assume the x and y position errors are known by measurement for the four 
positions. The position error for any point E inside the retangle ACBD can 
then be calculated by interpolation using 

E x = k-|C x + ( 1— k i ) Djj. 

Ey = k 2 Ay + ( 1“k 2 ) By 
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where the subscript denotes the component of the error and k-j , k 2 is the posi- 
tion of point E from C and A respectively. The interpolation is least 
accurate when E is farthest from the four reference points. This occurs when 
E is at the center or k-j = k 2 = 0.5. E at the center is the case illustrated 
by Figure 19. 


Analysis of the 157 points in Figure 19 showed that and 
Ey calculated by interpolation different from the measured values for E x and 
Ey by + 2 ym typically and + 4 ym at the greatest. 

Assume the interpolation error is proportional to the area of the 
ACBD cell. Then a + 0.5 ym interpolation error can be achieved using 
( 2)2 x 0.5/4 = 0.25 mm^, g.5 x 0.5 mm cell. The 12.5 x 12.5 mm will be 
mapped to a sufficient density using (12.5)2/(0.25) = 625 points. 

Changes in target light level and background do not directly 
affect the accuracy as long as a minimum signal-to-noise ratio exists. 

Section 3.5 shows SNR = 30 for +0.5 ym error. A signal-to-noise ratio test is 
included in the target acquisition routine to ensure sufficient target inten- 
sity for the desired measurement accuracy. The SNR test would add about 0.1 
second per target to the acquisition time. 

In summary the camera, lens and target intensity contributions to 
the loss of measurement precision are: 

Deflection Step Size: 

Temperature Coef of FOV Size: 

Temperature Coef of Center Of: 

Mechanical Package: 

Corrected Lens and Deflection 
Distortion: 

Signal-to-Noise Ratio: 


+0.4 ym 

+0.12 ym/°C at edge 
0 at center 

fset: +0.12ym/^C 

+0 . 5ym 


+0 . 5ym 
+0 . 5ym 
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3.6.2 Target dynamic motion effects upon system precision . ** The targets 

are moving while the cameras attempt to measure their positions. Section 3.3 
cited some of the amplitudes and frequencies that might be expected. Figure 
20 shows an envelope of possible amplitudes and frequencies for a variety of 
models. The amplitude is given in urn at the faceplate of a camera. 

A scale along the right side shows the number of tracker cycles needed to null 
a given amplitude offset to + 0.5 ym (Figure 18) and produce a position 
measurement. Diagonal lines labeled lOOg and 200g are lines of constant acce- 
leration at the model. They are calculated using g =• 4u 2 f 2 M x A/9800 , Mjj = 100. 
Their purpose is to provide some perspective on the vibration levels being 
used in the analysis. 

The tracker system measurements are time samples of the target 
position. The objective of the process is to obtain a sufficient number of 
samples to permit the user to estimate, in a statistical sense, the average 
target position. The model will be perturbed by random turbulence in the 
airstream. The turbulence will excite all of the dominant vibration modes of 
the model and modes of the model-sting combination. The vibrations will be 
damped and have dominant frequencies but their occurrence depends on the 
nature of the turbulent input. The time samples of position by the tracking 
system are assumed to be uncorrelated with the vibration. 

There is a period of time called the dwell time T. During this 
time the tunnel flow conditions and model angle of attack are held constant. 
According to NASA, a dwell time of 2 seconds is expected. This finite dwell 
time sets limits on the low frequency vibrations that the user can tolerate 
even through the tracker system can quite easily measure positions in the pre- 
sence of such frequencies. The lower limit is of the order of 1 Hz to 10 Hz. 

The time required to make a position measurement is 130ys. An 
additional 40ys elapses if the same target is repeated or 70ys if a new target 
is selected. The measurement is actually composed of two 45 ys submeasurements 
separated by a 40ys computation interval. For a target with vibrational 
motion in a direction 45° to the track scan axes, the x and y position values 
are taken at slightly different image locations. During the 85ys time inter- 
val after x is measured and until y is measured, the x position will have 
moved at most 

X = {2 A sin (2?ifAt) 

' o 


■f 

Ao 

X 

8 

100 ym 

0 . 3ym 

20 

iCOym 

0 .8ym 

50 

100pm 

1 .9ym 

200 

lOym 

0 . 3ym 


These errors go to zero as the motion becomes parallel to the scan axes , 
leaving the residual errors described for the camera and optics in Section 
3.6.1. 
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So far it has its been shown that the vibrational amplitudes are 
within the system capture range, that the measurement is made in a short time, 
I30y3 and that the various errors and uncertainties at the faceplate are of 
the order of 0.5 pm each. 

Several questions still remain. How many targets can the system 
handle, how does the precision depend on this number, how many measurements of 
a particular target are made, and what observing strategy should be used. 

When sampling a sine wave vibration of frequency f, two time 
intervals are important. The first one is t(n) * 130 ys (n) + 40 ys (n-1) 
where n *= 1. This is the time used to make a measurement of the vibration 
amplitude. While n cycles may be needed to converge to null, it is only the 
time involved in making the (last) position measurement that is important. In 
sampling systems this is called the aperture time. The second time interval 
of importance is T (N,n) ■ Nt(n) + (N) 70 ys. This is the time between 
samples where N targets are scanned and each target requires n cycles for 
convergence. 


The amplitude response of this sampling process is 


A » 


A 


o 


sin 


(TTft(l)) 

irft (1) 


(irfT(N,n)) 
sin IrTTWn) 


where 


A 0 is the peak-to-peak vibration amplitude, and 


f is the vibration frequency.. 

The first term, sin (u-ft ( 1 ) )/irft( 1 ) can be set equal to 1 since for all ampli- 
tudes and frequencies shown in Figure 20, its value is greater than 0.998. 

The argument TrfT(N,n) can be rewritten in terms of a normalized number of 
targets and in a manner showing the dependence upo.i vibration frequency and 
amplitude . 


Let 


irfT(N 0 ,n) = ir 

fT(N 0 n) = f N 0 (t(n) + 70 ys) « 1 

N = 1 1 

° f (tin) + 70 ys) 

= 1 1 

f (n(130 ys) + (n-1) (40 ys) + 70 ys) 

N 0 is a function of frequency and n, the number of cycles to converge to null 
+ 0.5 ym. The number of cycles n is , by Figure 18 related to the amplitude of 
the vibration. 
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The Figure 21 shows the plot of the normalized amplitude response 
of the system versus a normalized number of targets. The normalization factor 
for the number of targets is selected based upon the highest vibration fre- 
quency to be measured. Several values are calculated and listed in the top 
table of Figure 2 1 . 

If it is necessary to measure a broad vibration spectrum having 
large amplitudes and high frequencies, the number of targets will be 
restricted. The top table and graph of Figure 21 show if f - 200 Hz, Aq «* 10 
Un and n * 2 cycles to converge- a system using 5 targets will measure an S put 
peak-to-peak amplitude. If the sampling rate and number of targets is 
selected on the basis of 50 Hz as the maximum frequency, but Aq * 125 U ra, the 
system will measure 100 pm peak-to-peak and accommodate three targets. The 
chief reason for the reduced number of targets is the number of cycles needed 
to converge. More about that later. 

The cases cited have been extreme cases. They permit long term 
examination of relatively few targets . It is not always necessary to spend 
all 2 seconds of the dwell time looking at the same 3 or 5 targets. It is 
only necessary that no more than 3 to 5 high frequency or large amplitude 
targets be interleaved in the sampling process. At 2700 measurements per 
second, each of 5 targets has 54 samples in 0 . 1 second. The other 1.9 seconds 
could be devoted to other groups of 5 targets, provided the amplitudes of the 
low frequency vibrations, 1-10 Hz, are low. Should there be concern about the 
1-10 Hz content of the vibrations, consider 10 blocks of 5 targets, each block 
of targets is observed for 0.1 secon; a total of 1 second. The 10 blocks are 
then examined twice for a 2 second total. The point to be made is that the 
observing strategy depends upon the vibration spectrum. This aspect merits 
further examination. 

If the system is measuring 50 Hz, 125 ym vibrations at the calcu- 
lated 820 measurement per second rate, and if there should be a vibration pre- 
sent whose frequency exceeded 820/2 Hz, the RMS amplitude of this vibration 
will appear as a noise source added to the various other sources. This aspect 
of the target position sampling is similar to other sampled data systems. The 
input signal must be bandwidth limited or else aliasing of the high frequency 
components occurs which addes to the in-band noise . 

The system presented so far has not been completely optimized. 
Operating conditions have been representative of fast, general purpose image 
dissector cameras . The tracker algorithm is a very simple one . In Section 
3.5.1 it was mentioned that the error detector curve was non-linear. It is 
this non-linearity that costs most of the system speed at large vibrtion 
amplitudes. The tracker is slow to converge to null; seven samples are needed 
in the case of 125 y m offset. 

More light at each target allows the linearization of the error 
detector curve by reduction of slope through null. A 2X reduction in the 
number cycles needed to converge various offsets seems possible. The lower 
half of the Figure 21 table is based on such a 2X reduction 
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3.6.3 Camera dynamic m otion effects upon system precision . - For a 

single camera and small '•ude vibrations, the changes in apparent target 

position due to camera vibrt-... ' cannot be distinguished from changes due to 
the vibration of the target. 


Consider a target at Point P-|. The vector between the image at 
the camera and P-j pivots at the focal point. 


Measurement of P-j in the presence of small amplitude vibrations ‘ 
about p.j is the purpose of the system. (Small, in this case, is in comparison 
to the distance between F and P i . ) 

If the target moves a small distance AP perpendicular to FP-j , the 
corresponding change at the image plane is 

Ai = AP 

b 

where a and b are the magnitudes of I-jF and FPi» 



If the camera moves in a plane perpendicular to I-jF by an amount 
f the focal point F moves . The corresponding Ai change is 




Equal magnitudes of uncorrelated target or camera motion lead to 
the same image motion. Therefore target motion plus camera motion must be 
vector summed and the size of the required capture range determined from the 


sum. 


Magnitudes and frequencies of camera vibrations were unknown at 
the time of this study. There will however be vibrations. 

3.6.4 Seeing conditions . - The light from a target passes through the 

air flow in the tunnel, an optical window in the wall of the test section, and 
the window on the pressure isolation container for the camera. All three 
potentially effect the image quality and the measurement process. 

The effects of the windows can be minimized if they are thin flat 
and not allowed to move relative to the optical axis or camera to model vec- 
tor. Motions in the plane of the windows have no effect since there is no 
change in the optical path length or angles of the rays . 

Tilts of the windows can be caused by pressure changes, or tem- 
perature changes . Tilts of the cameras relative to the window have a similar 
effect. 
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For a ray passing through a 10mm thick quartz window at an angle 
of 35° to the normal, an angular change of +0.5 degree will cause an apparent 
+65 pm shift in target position at the model. The Figure 22 illustrates this 
point. Note that the apparent target shift is proportional to the window 
thickness . 


The air in the tunnel also acts as an optical element. Even 
though the index of refraction is very close to 1.0, the interface at the win- 
dow has an effect. The wind tunnel static pressure can change from 1 to 9 
atmospheres and the temperature from -195°C to 95°C. The index of refraction 
n^ of a gas varies as 


n 

t 


(n -1) 
o 

(1-0. 00366t) 


P 

760 


where n Q = 1.000298 for nitrogen at 0°C and 76:im-.n Hg. 


as 


For the temperature and pressure extremes, the n^ value calculates 


n t ( -195°C, 9 atm) = 1.00937 

n t (+ 95°C, latm) = 1.000221 

For rays passing at 35° through the 10mm thick window in the test 
section wall, the change is 47ym in apparent target position at the model. 
Again this shift is proportional to window thickness. 

There is a variation in this window induced error over the field 
of view since the angle is 35° only at the center. The view changes from 
18° to 52° at the edges of the field of view. 

The effects of the temperature and pressure gradients in the vici- 
nity of the test section window and near the model have not been analyzed. 

3 . 7 Use of Tracker Cameras for TV Viewing of the Tunnel Interior 

Image dissectors were originally conceived and developed for tele- 
vision camera tube use. They have been largely replaced by camera tubes 
having image storage (vidicons, etc) because of the greater light efficiency 
of frame storage devices . 

Image dissectors are mainly used where slow scan is needed; as in 
narrowband width image transmission systems. The image dissector signal 
current is independent of the scan rate. The signal current in storage camera 
tubes is proportional to the scan rate. 

For television purposes, the image dissector camera resolution is 
set by the aperture in the tube. The nominal size being considered for the 
tracker system is 200ym x 200ym square within a field of view of 10 x 10mm 
square. This results in a TV picture with 100 percent modulation at 50TVL/ 
picture height and zero modulation at lOOTVL/picture height; a rather low 
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resolution image. The typical closed circuit TV resolution of 250TVL/picture 
height at 50 percent modulation can be achieved with an aperture 50 x 50pm 
square. Based upon the analysis of the tracker system, it is doubtful that 
the tracker would work using a 50Mm aperture; the capture range would be 
severely reduced. 

Multi-aperture image dissector tubes have been built in the 25 and 
44mm sized tubes. There is no basic reason that a dual aperture version of 
the 17mm tube could not be built. If one of the apertures is 50 x 50pm while 
the other is 200 x 200pm, the needs of the TV and tracker modes could both be 
satisfied. However, it can be shown by the following calculation that the TV 
mode must be slow scan, minutes per frame; or else require high tunnel illumi- 
nation, 1000ft cd. 

The highlight signal-to-noise ratio of a TV picture is required to 
be at least 35:1 BMS . At this level the image will begin to show noise 
breakup. The signal-to-noise ratio of an image dissector camera in the TV 
mode is proportional to the square root of the number of photoelectrons 
passing through the aperture in a time t = 1/2BW; BW = video bandwidth. 

e K SNR 2 

Ipixel = — — amperes 

e = 1.602 x 10“1® Coulomb/electron 


K = 2.5 noise factor for image dissectors 


Ipixel = a,j°2g-.19M.2.5) .(352) 

The time t is related to the number of pixels in a TV frame and the frame 
time. Let the 10 x 10mm field of view be scanned by the 50pm aperture in 25pm 
steps. There are then 1.6 x 10^ pixels per field of view. Furthermore, let 
there be a 5 percent allowance for line and frame retrace . The frame time T 
is therefore 


or 


(1.60 x 10 5 )At 
1 " (0.95)2 

At = (0. 95) 2 T 

1.60 x 10i» 


The pixel current is therefore 


Ipixel = 


(1.602 X 10~ 19 ) (2.5) (35) 2 (1.60 x 10 5 ) 
(0.95) 2 T 

87 x 10~ 12 

amperes 


Image dissectors are rated in terms of maximum total photocurrent. The maxi- 
mum values is 25pA for the 17mm tube. 
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The total pliotocurrent is 


Ipc = (1.6 x 10^) Ipixel 

- 13.9UA < 25 yA 

T 

Therefore, the shortest frame time is about 0.5 second or 2 frames/second. 

This TV picture will have noticeable flicker for the viewer. 

In order to obtain 25yA of photocurrent, the typical 200 A/lumen 
photocathode needs a light level of 0.125 lumen or 115 ft cd. This is bright. 
Scene illumination is 4f #2 greater or 3600 ft cd for f/2.8. 

The conclusion is : use a vidicon . 

3 . S Demonstration 

A single camera operating under control of a 16-bit mini-computer 
was used to demonstrate multi-target position measurement and tracking. The 
demonstration showed the ability to search and acquire five targets, to track 
the movement of these five targets while they were moved at a rate and in a 
manner simulating a change in the angle of attack of a wind tunnel model, and 
to measure the position of the five targets under conditions simulating the 
target vibrations described in Section 3.3. 

Since the demonstration used laboratory equipment, a general pur- 
pose camera, and a tracker function implemented in software, some scaling of 
the conditions were necessary. 

3.8.1 Field of view . - The demonstration field of view was 4 x 4mm com- 
pared to the desired 10 x 10mm. The smaller field of view was selected in 
order to achieve 0.5ym deflection increment using the available 13 bit D/A 
converters. Use of a 15 bit converter permits the entire 10 x 10mm field of 
view to be utilized at 0.4ym deflection increment. 

3.8.2 Targets . - The five point targets were pinhole apertures in a 
metal plate. The targets were back illuminated by a tungsten light source. 

The image of the test pattern was demagnified onto the faceplate of the image 
dissector camera using a 55mm f/2.8 lens. The background brightness was mini- 
mized to avoid the added task of writing a background suppression routine in 
the target acquisition software. (See Section 3.2.3.) 

The test pattern was mounted on a x/y table that permitted it to 
be moved about simulating angle of attack changes . 

3.8.3 Target acquisition . - A program subroutine caused the camera to 
begin raster scan of the field of view in 50ym steps. If the video signal was 
above an absolute threshold, the target list was checked for the current coor- 
dinate or one within a specific range. If no target was listed, the present 
one was added. If the target was listed, the scan skipped a safe distance 
beyond the target and continued until the raster was completed. The time to 
acquire the five targets was about 5 seconds. 


Five 25 pm Apertures at Locations Indicated 
Used with 4:1 Demagnifying Lens 





DIMENSIONS IN mm 


Figure 23. Demonstration Test Pattern 


3.8.4 Track and measurement . Upon completion of the acquisition, the 

program jumped to the track routine. This routine measured each target posi- 
tion in sequence, updating the target list. The routine was adaptive. It 
would remain at a target, counting cycles, until the process converged to an 
operator specified limit or until five cycles occurred. Upon convergence or 
five cycles the routine moved to the next target. 

The adaptive ability allowed the tracker to follow very rapid 
linear motions of the targets; velocities much greater than would occur during 
an angle of attack change. This mode assures that no targets would be lost 
while translating the test pattern rapidly about the field of view. 

The tracker action was a time scaled version of the timing diagram 
shown in Figure 8 . Figure 24 shows the comparison between the demonstration 
of the tracking algorithm. All software was written in FORTRAN to ease the 
software development task. There is also a sample called display. This was 
used to show, by means of a 'scope in parallel with the camera, the converged 
estimate of the target position. 

Four time scales are involved . They are shown in Table 2 . The 
five demonstration targets are scanned in a period equal to about 160 targets 
in the design concept system; or 33:1 ratio. 

3.8.5 Target dynamic simulation . - A subroutine was used to generate a 
table of vibrating target position versus time coordinates . The amplitudes 
and frequencies used are given in Table 3. The frequencies were scaled to 
present the same time scaled motion to the demonstration system as the actual 
frequencies would present to a real time system. 

The table of vibrating coordinates was then used to disturb the 
target coordinates contained in the target list before the target coordinates 
were sent to the tracker. This creates the illusion that the target has moved 
since last examined. The tracker then must reconverge to the null position. 

Of the conditions listed in Table 3, the tracker occasionally 
experienced difficulty with the third one, the yaw vibration of 126ym x lOym 
at 50 Hz. Using the 125x rule, this corresponds to 15mm x 1.3mm at the model. 
It was not clear if this was a dynamic limit or some unresolved problem with 
the demonstration hardware. The 126ym is inside the capture range. 

Table 2. Comparison of Stereo Slectro-Opical Tracker System 
Time Scale with the Demonstration System Time Scale 



Demonstration 

Stereo E-0 


System 

System 

Sample Time 

200ys 

20ys 

Cycle Time 

2600ys 

130ys 

Intercycle Time 

2900ys 

40 s/70ys 

TOTAL for 3 

18 . 1ms 

540ys 

cycles at 



1 target 
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and Laboratory Demonstration System 


Table 3, Vibration Simulation Frequencies and Amplitudes 


Extremes are 


Pitch: 

13,4PM 

3,0pm 

50 Hz 

Roll: 

25,7pm 

7,7pm 

50 Hz 

Yaw: 

126,0pm 

10,0pm 

50 Hz 

Wind Dynamic: 

78,0pm 

38,0pm 

20 Hz 

Wind Dynamic: 

7,8pm 

3,8pm 

200 Hz 


All Dimensions are Peak-to-Peak at the Photocathode. 


4.0 


CONCLUSIONS AND RECOMMENDATIONS 


This study and report provides an analytical and experimental 
basis to support the concept of a stereo electro-optical tracker system to 
measure multiple targets attached to aircraft models in the wind tunnel at the 
National Transonic Facility. The study showed given 0.8m x 1.0m field of view 
and 150 um precision at the model/ 5-50 point targets can be acquired and 
tracked. Furthermore, the intensities of the targets, tunnel background illu- 
mination, and system data rates are within reasonable bounds. A demonstration 
showed a time scaled simulation of a five target tracker operating at the 
expected electro-optical signal-to-noise ratio and desired precision. The 
analysis and demonstration also showed that targets vibrating at frequencies 
of 8, 50, and 200 Hz can be located and measured. There is a trade among the 
number of targets per second, the vibration amplitudes, target brightness and 
contrast that must be considered. 

Several important areas remain to be examined or specified. It is 
these areas that introduce an element of risk to the development and use of 
this system concept. 

Targets and their attachment to the model are critical to the suc- 
cess of this concept. This study proposes that LEDs embedded in the model be 
used. An alternative is fiber optic light pipes, ends polished flush with the 
model . 


The effects of the distortions in the camera and camera optics on 
accuracy needs examination from the viewpoint of the stereo-optic transfor- 
mation. These distortions can be corrected? but knowledge of the required 
degree of correction for all sources is needed,. 

The seeing conditions between the target and the model are impor- 
tant. Do the shock waves, turbulence and other pressure gratients degrade the 
accuracy of the system? 

System calibration needs further examination. The study has con- 
sidered only some of the basic aspects and the general approach. 

Lastly, it will be important to know what kind of a thermal and 
vibrational environment the cameras will experience inside the tunnel pressure 
vessel. This need is not just for questions of mechanical packaging. 
Vibrations will contaminate the measurement process . The study shows the 
system can measure vibrating targets on a vibrating model. System speed and 
number of targets depends strongly on the vibration amplitudes and frequen- 
cies. What are the limits on camera-to-camera and model-to-camera motion. 

How does the user/experimentor extract useful information about true target 
position from the three dimensional vibrationally perturbed data recorded by 
the system. 
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4.1 


Credits 


The writer has found this concept of a stereo electro-optical 
tracker system a most facinating one to examine. The help of David Gray of 
NASA Langley was most appreciated. He answered many questions about the tun- 
nel environment and didn't mind repeated requests for vibration data. The 
software for the demonstration was written and debugged in about three weeks 
time by Doug Renselle of SCS, Inc. His contributions greatly aided the 
demonstration. 
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APPENDIX 


The Image Dissector as an Optical Tracker 
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THE IMAGE DISSECTOR AS AIM OPTICAL TRACKER* 


APPLICATION NOTE E16 


1.0 INTRODUCTION 

Although the image dissector was probably the first all-electronic optical scanning device to be used (Ref. 1) 
for generating a video television signal, its first reported use (Ref. 2) as an optical tracking device did not occur until 
1960, Since that date it has enjoyed a steadily increasing popularity as a specialized optical scanning device, especially 
for tracking star (Ref. 3) and laser (Ref. 4, 10) images. 

2.0 BASIC PRINCIPLES 

To understand why the image dissector is so well adapted for optical image tracking, consider first its 
mechanical analogy shown in Figure la. In this mechanical analogy an opaque mask with an aperture is moved in the 
plane of an image, until the desired image flux falls through the aperture and onto a suitable detector, such as a 
photomultiplier tube. By determining the precise mask position at which this aperture penetration occurs it is then 
possible to determine the position coordinates of the image. 

More importantly, perhaps, from the standpoint of image dissector operation, once the image is located, its 
subsequent motion can be tracked by starting a series of small incremental motions of the aperture, without the 
necessity of continuously, and inefficiently, scanning the total viewed field. It is this ability of the image dissector 
to adapt its own scanning motion to that of the image to be tracked, which distinguishes it from the majority of other 
optical trackers. 

The electronic image dissector shown in Figure lb can perform exactly the same function as the mechanical 
dissector, only much faster and without backlash. In this case, the optical image is converted first to an electron 
image, by a suitable combination of photocathode and electron lens. The electron image is then moved, in a nearly 
Inertialess manner, over a fixed mask containing an aperture. Behind the aperture is a high gain electron multiplier 
which converts each entering electron into a large, often individually detectable (Ref. 5) burst of charge in the 
dissector anode circuit. 

3.0 TYPICAL IMAGE DISSECTOR DESIGNS 

Commercially available image dissectors vary primarily in the methods selected to focus and deflect the 
electron image, as well as in the diameter of the input optical image which can be accommodated. Figure 2 shows the 
basic configuration used in three popular types of image dissectors. 

* Paper presented by E. H. Eberhardt at the Optical Tracking Seminar of the Society of Photo-Optical In- 
strumentation Engineers, El Paso, Texas, January 18-20, 1971. 
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At the top, Figure 2a, is a Vidissector'fi' (Footnote 1), image dissector, an advanced version of the earliest, 
magnetically focused and deflected dissectors. A fine pre-acceleration mesh is used in close proximity to the photo- 
cathodc, to accelerate the photoelectrons up to their full velocity as quickly as possible, followed by an electric-field- 
free drift and deflection space prior to the defining aperture plate. This tube is designed to operate satisfactorily in the 
same type of magnetic focus and deflection qoils used for vidicons and image orthicons. Such an image dissector was 
used as a star tracker in the first X-15 experimental aircraft (Ref. 6). Its properties as a star tracker have been 
described by Ostroff and Romanczyk (Ref. 7). 

In Figure 2b is shown the basic configuration of a hybrid image dissector with electrostatic focusing and 
magnetic deflection. This type of dissector, often called a “star tracker” tube, was used in the first reported dis- 
sector tracker (Ref. 2) and in other, more recent trackers (Ref. 8, 9, 10), including the star trackers used in Aerobee 
rockets and on the Lunar Orbiter spacecraft (Ref. 3). 

In Figure 2c is shown the basic configuration of a Reconotron all-electrostatic image dissector which was used 
as an azimuth star tracker on the Mariner series of spacecraft (Ref. 1 1). 

4.0 PERFORMANCE CHARACTERISTICS OF IMAGE DISSECTORS 

4.1 Rise for Qccultation) Distance 

When a small, “point-source” optical image is focussed onto an image dissector, as indicated in Figure 3, and a 
linear deflection signal is applied, the aperture can be made to move (Footnote 2) across this point source image, 
generating a pedestal-like output signal. Ideally this signal should have a sharp rise and decay characteristic, as in- 
dicated by the dashed lines in Figure 3 as the aperture occults the point source image, in practice, the rise and decay 
characteristics are somewhat rounded, with the degree of rounding conveniently describing the capabilities of a 
particular image dissector. Quantitatively, the rounding can be measured in terms of the aperture displacement, w, 
measured at the photocathode, necessary to change the output signal between the 10% and 90% response points, as 
indicated. The magnitude of the resulting rise distance, or occultation distance, w, is an important operational 
characteristic of an image dissector. 

Figure 4 shows how a finite rise distance is generated internally within an image dissector. The group of 
photoelectrons emitted from a point image on the photocathode, exit from the photocathode with a finite spread of 
energies (commonly averaging about 0.5 eV) and directed at various angles (commonly with an approximately 
Lambertian distribution). Thus these photoelectrons tend to spread out as they leave the photocathode, progressively 
destroying the image quality. This tendency to spread is counteracted by the focussing action of the electron lens 
which brings the bundle of emitted electrons back to a near-point-source size at its focal plane. It can be shown that 
the width-at-half-maximum of the resulting, approximately Gaussian, beam density distribution at the image plane is 
approximately equal to the 10%-90% rise distance, w. Thus, there are four essentially synonomous terms for the 
parameter, w, at focus, namely: 

w = beam width = beam diameter = rise distance = occultation distance 

Figure 5 shows what the total electron beam profile actually looks like in Vidissector image dissectors, for a 
point source input, and three possible magnetic focus conditions, namely: one loop, two loops and three loops of 
focus. Focus at the aperture can be assured by choosing a combination of the solenoidal magnetic field strength, 

I: Trademark, International Telephone & Telegraph Corp, 

2: Although the electron image actually scans over a fixed aperture in the image dissector, it is often helpful in 

describing dissector behavior to talk in terms of the mechanical analogy, where the aperture is moved over the 
optical image. 
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B n in webers/meter2; the photocathode-to-mesh voltage, V, in volts; the photocathode-to-aperture spacing, L, in 
meters; the photocathode-to-mesh spacing, L m , in meters; and the charge to mass radio, e/m, of the electron, in 
coulombs/kg, such that the following equation is satisfied: 


Bn 


27T T) 

L + L m 


2 mV !/ 2 


( 1 ) 


Insofar as both the electric and magnetic fields are uniform, and the aberrations introduced by the non- 
uniform deflection fields can be ignored, the beam width, w, at focus is given by: 

(V ro V zo ) V 2 

W a 1.23 L m - (2) 

where V ro is the average radial emission energy component of the photoelectrons, and V zo is the average axial 
component, in electron volts. 

For a typical Vidissector image dissector, such as an F4052, with L = 22.5 cm, L m = 0.32 cm, V =600 volts 
and 2 V ro s V zo s 0.3 eV (white light input to a multialkali photocathode), Equation 2 yields a value of 1.4 
micrometers for the beam width, w. This is approximately 0.003% of the total image field diameter which can be 
scanned in such a tube (4.2 cm), illustrating the excellent limiting “resolution” characteristics which can be 
achieved. Measured rise distances as short at 2-5 micrometers, have actually been observed in an F4052 image 
dissector camera. 


The maximum (off focus) beam diameter, D n , and the beam convergence angle, 0 (useful in estimating depth 
of focus) are given in Vidissector image dissectors by: 


nD n 3 0.53 (L + L m ) (V r0 /V) V 2 

(3) 

tan (0/2) a 0.83 (V r0 /V) l/ 2 

(4) 

Possible Tracking Modes 



One type of simplified “rosette” scan which can be used (Ref. 2, 3) for tracking a star or laser image with a dis- 
sector is shown in Figure 6. Here one applies successively, a small positive horizontal aperture displacement and 
return, a small positive vertical aperture displace and return, etc., as shown, with the star image initially known to be 
located somewhere within the aperture area. If the star image is centered within the rest position of the aperture, the 
output signal, as illustrated, will be a uniform symmetrical series of pulses, of rise distance, w. 

On the other hand, if the star is displaced, the output pulses will become assymmetrical in width, generating a 
fundamental component of the scan “frequency”, f 0 , whose amplitude is linearly proportional to the offset dis- 
placement of the image, and whose phase determines the direction of the displacement (see Figure 6). This type of 
simple, but efficient star tracker, combined with mechanical gimbdling to keep the stellar image within the viewed 
aperture area has been used (Ref. 3) to track stars within less than 1 micrometer displacement accuracy at the 
photocathode. For a particular optical system used, th ,s was equivalent to better than 10 arc seconds tracking ac- 
curacy. Note that the tracking information, to a first approximation is independent of the magnitude of the rise 
distance, w, and of the magnitude of the stellar flux, F, (assuming that the magnitude of the f 0 frequency component 
is referenced to the magnitude of the 2f 0 component). If so desired, the electronically generated image offset in- 
formation can be converted to a “DC” centering deflection displacement of the aperture. The aperture can then be 
made to electronically follow the motion of the star image within the field of view. 
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Four other possible tracking scan modes are shown in Figure 7. The first, Figure 7a, illustrates the azimuth- 
only type of tracker used in the Mariner spacecraft (Ref, 1 1). This tracking mode takes advantage of the fact that a 
long narrow, slit-shaped aperture can be built into the image dissector tube, with line scan permitting azimuth 
position information to be extracted without regard to the elevation of the star image. 

Figure 7b is a similar tracker, using an L-shaped aperture. Here, both azimuth and elevation information for a 
star image (or images) can be generated by use of a corresponding L-shaped scan. 

Figure 7c illustrates a somewhat more complex type of image dissector tracker. In this case a small aperture 
is rotated, or wobbled, with a small amplitude, such that image irradiance gradient information is generated. The 
center of the scan can then be made to move electronically at right angles to this gradient. The dissector then auto- 
matically traces out equal irradiance level contours. Such a scheme has been used in an eyeball tracker (Ref. 12) for 
tracking the pupil of the eye. 

And finally, a rather sophisticated image dissector tracker can be constructed whose internal aperture can be 
adjusted electronically to a specific shape, such as the letter, “A”, illustrated in Figure 6d. When such a custom- 
shaped aperture is then scanned over an image field, it will generate a large, sharply defined, “correlation spike” 
whenever it crosses precisely over an image area of the same shape. This sharp correlation spike can be used (Ref. 13, 
14) with great effectiveness to track complex image shapes against a non-uniform background. 

If so desired, dissector tubes may be supplied with multiple apertures, each with its separate read-out electron 
multiplier, for special optical tracking applications (Ref. 15). 

4.3 Signal-to-Noise Ratio 

Usually, only one type of random noise is encountered in image dissector trackers, namely the inherent 
statistical fluctuations of the photoelectrons emitted from the photocathode. Dark noise, i.e., random residual 
noise generated by the tube itself, and by the external amplifier circuits, in the absence of input flux, is nearly 
always negligible. This is a direct benefit of the large (> 1(F), nearly-noise-free charge amplification present in the 
electron multiplier of the image dissector. 

To calculate the expected magnitude of the photoelectron noise and show how it limits the ability to find 
and track images, it is convenient, and especially instructive of image dissector behavior, to assume that the dissector 
scans by making a series of discrete image intensity samples. In this sampling process, illustrated in Figure S, the scan 
stops on each aperture-sized element for a short time, A t, called the dwell time. It can then make a “count” 
(Footnote 3) of the total number of photoelectrons collected by the aperture. 

When only a uniform background irradiance is incident on the photocathode, the resultant output count, n^ 
will be given by 

nj, = LaSjj E A t (5) 


3: While the image dissector can actually count individual output anode pulses corresponding to individual photo- 

electrons entering the aperture, it is more usual in practice to accumulate these anode charge pulses and generate an 
average output current. 
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where 


L = background irradiance, incident on the photocathode in (flux units) m'^ (Footnote 4) 

a = aperture area, in (Footnote 5) 

S\j = photocathode responsivity ratio, for the background irradiance flux, in electrons ■ sec’* ■ (flux unit)'^ 

(Footnotes 4 & 6) 

E = counting efficiency of the electron multiplier, in output counts per emitted photoelectron collected 

by the aperture area a. 

At = dwell time, in seconds 

Of course, if this background-generated count, n b , were exactly constant each time it was made it would not 
interfere, in any way, with the subsequent ability of the image dissector to And image areas of greater or lesser input 
flux. But it is not, and cannot be constant because of the inherent random emission of photoelectrons. It will there- 
fore have the usual Gaussian variance, or noise content, given by: 

n b l / 2 = (LaSbEAt) 1 / 2 (6) 

The noise problem in image dissector trackers under these conditions then resolves itself into whether or not an 
observed change in the count for a particular image element is an accidental change, according to this expected 
noise variance, n b V 2 , or whether it is due to true change in the average flux input. For a point image input (i.e., a 
star or laser image) the increased count, s, due to the incident stellar flux, F, is given by: 

s = signal count = F Sf E A t (7) 

where F = stellar flux incident on the photocathode, in flux units (Footnote 4) Sf = photocathode responsivity 
ratio for the stellar flux input, in electrons • sec*l • (flux unit)'l (Footnotes 4 & 6). 

Note that the signal count, for this special case where the signal flux is confined to an area smaller than the 
aperture, is independent of the aperture area, a. 


4: The units used to measure flux are optional. Photons/second, watts, or lumens, are the three flux units most 

commonly selected. 

5: In order to automatically correct for the sometimes-unknown electron optical demagnification of the electron 

lens, aperture dimensions are nearly always expressed in terms of their effective magnitude measured at the photo- 
cathode surface. 

6: The magnitude of the photocathode responsivity ratio, S, in electrons > sec*! (incident flux unit)-l, depends 

on the flux units selected and on the spectral distribution of the flux input. Methods of computing this ratio for 
various flux units and various spectral distributions have been described and tabulated (Ref. 16). 
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An effective signal-to-noise output power (Footnote 7) radio SNPRtrack for image dissector star trackers can 
therefore be written down directly as: 


F 2 Sf 2 E 2 At 2 

SNPRtrack » " (F 2 /L) (Sf2/S b ) E (At/a) (8) 

L aS b E A t 

The importance of maintaining high stellar flux, F, high cathode responsivity to this stellar flux, Sf, and high 
electron multiplier counting efficiency, E, is evident from this equation. 

However, changes in the dwell time A t, and aperture area, a, must be interpreted with care, since there is 
usually an interdependence between the aperture area, a, and the dwell time A t. This can be illustrated for the 
special cose where a square aperture image dissector scans a total field area, A, in a total search time, T, with no 
scan overlap and negligible fly back time. In this special case the ratio, At/a equals T/A, and the signal-to-noise ratio 
becomes: 

SNPRtrack = (F 2 /L) (Sf 2 /S b ) E (T/A) (8a) 

In this alternative form, and for this scan pattern, it can be seen that the signal-to-noise ratio is independent of 
the selected aperture area. Thus in certain star tracking modes one is free to choose any desired aperture size, and 
thus to locate the star image as precisely as needed, without loss of signal-to-noise ratio. 

This rather surprising deduction does not hold when the star image becomes comparable to the aperture size, 
or when the signal is derived from changes in the input fleet density level, as in TV-type scanning, where the signal is 
linearly proportional to aperture area, a. For these latter situations, the more conventional signal-to-noise ratio re- 
lationships (Ref. 17) will apply. 

If a current measuring external circuit with a noise bandwidth, Af, is used instead of a counting circuit, the 
dwell time, AT, is usually replaced by: 

At a (2 A f)’ 1 (9) 

and the counting efficiency, E, by: 

E a k' 1 (10) 

where k is the noise factor (Ref. 18) of the electron multiplier. Counting efficiencies of 0.85 to 0,95 output counts 
per emitted photoelectron are easily achieved in image dissectors with well-designed first dynode geometry (Ref. 5). 
This is reduced by the electron transmission ratio of the mesh (typically 0.6 to 0.65) if a mesh is used, as in Vidissector 
image dissectors, between the photocathode and the first dynode. 


7: Image dissector signal-to-noise ratios are also often written in terms of the ratio of output signal current to 

output noise current. 
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5.0 


ADVANTAGES. DISADVANTAGES. AND PROBLEMS OF IMAGE DISSECTORS 


The optical tracking systems designer, faced with the problem of selecting an appropriate tracking method, 
needs to know both the advantages and the disadvantages of the image dissector, as well as some of the problems 
which can arise. Unfortunately most of the important problem areas concern specific tube types, and often specific 
tube samples, so it is difficult to make generalized conclusions. Nevertheless some useful statements can be made. 

5.1 Advantages of Image Dissectors 

5.1.1 Fast Deflection: The low inertia of electrons allows fast deflection. Limits are more often established by 
the external deflection circuits, 1 M sec jump times between adjacent image elements and 100 n sec jump times 
across a full image diameter have been achieved, and improvements are in progress. 

5.1.2 Adaptable Scan: The image dissector can easily be switched on command, from one scan mode to 
another, e.g. from search scan to track scan or from search scan to contour-tracing scan, or from jumped scan to 
continuous scan. 

5.1 .3 High Tracking Mode Sensitivity . When the aperture of an image dissector is actually observing an optical 
image to be tracked, it enjoys very high effective sensitivity, equal to or better than the best photomultiplier tubes. 
No difficulty is experienced in counting individual photoelectrons, and the spurious dark count encountered is ex- 
tremely small (typically a few counts/second or less), 

5.1.4 Linear Input-Output : The output from an image dissector is linearly proportional to the input flax 
over many orders of magnitude. Thus the image dissector can be used as a precision microphotometer simultaneously 
with scan tracking (Ref. 19). This property is also often useful in deriving optical intensity information for a com- 
puter input (Ref. 20). 

5.1.5 Wide Bandwidth : The image dissector can follow intensity modulation of the input flux seen by the 
aperture from DC up to many megahertz. With special electron multipliers this can be extended into the gigahertz 
region, making the image dissector especially attractive as a dual mode sensor for optical tracking and wide band- 
width optical communications systems (Ref. 4, 28). 

5.1.6 No Gear or Rearing Backlash: Since the aperture as well as the deflection and components can be per- 
manently attached to the tube envelope, there are minimal opportunities for boresight offset errors to occur. 

5.1.7 Rugged Construction: The inherent simplicity of image dissector construction permits the tubes to be 
ruggedized. No difficulty has been experienced in meeting the mechanical requirements for rocket and missile 
launch; in fact at least one image dissector tracker tube has survived an Aerobee rocket crash landing and was sub- 
sequently re-used. 

5.1.8 Long Life: Dissectors have now been in essentially continuous service for over 30 years. The absence of a 
short-lived thermionic cathode, used in other TV camera tubes, make this performance possible. 

5.1.9 Damage Resistance: The dissector is inherently resistent (Ref. 21) to many types of damage such as 
momentary flux input overloading, nuclear radiation exposure, accidental over-voltaging, etc. For example, a 
Vidissector image dissector in the ATS-F spin stabilized satellite (Ref. 22), has withstood periodic direct solar ex- 
posure through a high speed optical system, with no sun-shuttering and with no removal of voltages, for over three 
years in space without measurable degradation. 

5.1.10 Fast Turn-On : An image dissector is ready to go within a few milliseconds, at most, after application of 
the proper operating potentials. 
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5.1.11 Wide Spectral Response: Photocathodes for image dissectors can be formed with spectral responses 
ranging from the far ultraviolet region, to the near infrared ( 1 J Q nm to 1200 nm), 

5.2 Disadvantages and Problem Areas 

. %•> » *»< H** 

5.2.1 Low Search-Mode Sensitivity : Since the image dissector wastes alf image photoelectron information 
which does not happen to pass through its aperture at any instant of time, it is inherently less sensitive than “storage” 
type scanning tubes such as the image orthicon and vidicon, for the wide field search mode, where a large number of 
image elements must be searched before the image to be tracked is located. 

5.2.2 Non-uniform Response (Shading, Blemishes, etc.): The input-output responsivity of an image dissector 
wili generally vary, at least by a small amount, from one point to another over the viewed field. If the change is 
gradual, it is called “shading", with variations within + 10% or so over the image field being typical. More serious, in 
many applications, are the rather sharp variations called blemishes which can occur within one or more picture ele- 
ments. The number and amplitude of these blemishes is often an individual tube-sample characteristic. However, as 
in other TV camera tubes, blemish numbers can be kept small by proper quality control during tube manufacture. 

5.2.3 Off-Axis Defocussing : When a peripheral image element is deflected into the aperture of a dissector, the 
effective aperture size and shape, and the rise distance, w, may change. This can be, and commonly is correctable in 
practice by applying a dynamic focus signal to the tube focus circuits in synchronism with the deflection. Such 
dynamic focus signals are usually easy to appl> . The residual defocussing, after correction, is partially limited by 
tube design and partially by the external focus and deflection coil design. In Vidissector image dissectors, with cus- 
tom wound coils, residual rise distances as low as 3-5 «m after dynamic focusing, have been achieved over the full 
usable input image area. 

5.2.4 Off-Axis Distortion : The actual deflection produced by applying a deflection signal is not necessarily 
identically proportional to the deflection signal. The resulting image distortion is especially serious in electro- 
statically focused tubes, which are really satisfactory only for near-axis tracking. In magnetically focused Vidissector 
tubes the image distortion is largely a function of the coil design geometry, with aperture displacements within less 
than 1% of their linearly predicted position being readily achieved. 

5.2.5 Boresight Hysteresis: Although the image dissector has no gear and bearing backlash problems, there 
can be an equivalent problem introduced by hysteresis in the deflection circuits (failure to reproduce exactly the 
same aperture position for a given applied deflection signal). In magnetically focused tubes this can be caused by 
magnetic hysteresis in the magnetic core material or magnetic shielding material. In electrostatically focused tubes it 
can be caused by electrostatic charging of internal insulator surfaces. With proper tube and coil design these effects 
can often be reduced to negligible values. 

5.2.6 Scan Delay : There can be a finite delay between the time a deflection signal is applied and the time the 
aperture reaches its new position. This can be due to eddy current losses occurring in the metal parts in magnetically 
deflected tubes, or to finite insulator charging time in electrostatically deflected tubes. This problem can be 
especially troublesome in random jumped scanning (not usually encountered in optical trackers) in which the 
aperture must move very rapidly over large distances and then be suddenly stopped. 

5.2.7 Photocathode Fatigue : With continuous intensive current density loading of a photocathode, it is 
possible to gradually reduce its responsivity, especially in the red region of the spectrum. By proper tube design, this 
photocathode fatigue can be greatly reduced, with total photocathode emission charge densities ranging up to 
several coulombs/cm2. Barring tube breakage, sudden, catastrophic failure is almost never encountered. 
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5.2.8 Gain Stability : The high gain properties of electron multipliers make them inherently sensitive to the 
magnitude of the applied voltage and to small chemical and physical changes in the dynode surfaces. Thus the 
applied voltage should be regulated, and the temperature held within reasonable limits, such as + 10°C, to achieve 
stable gain. 


6.0 SUMMARY 

It con be seen that the image dissector is nearly an ideal device for certain optical tracking applications. It can, 
and does, compete with mechanical trackers in some areas, and in other areas, with more conventional 
raster-scan-only television camera tubes. 
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